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^ Based on observations obtained at the ESO VLT 
on Paranal, Chile, under programmes 63.1-0015 
and 65.1-0135. 

3 As defined in Stolte et al. 2004 (Paper I), 
we refer to the central cluster in NGC 3603 as 
"NGC3603YC" to avoid confusion with the ex- 
tended HII region. 
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ABSTRACT 

Based on deep VLT/ISAAC JHK photometry, we have derived the present- 
day mass function of the central starburst cluster NGC 3603 YC (Young Cluster)^ 
in the giant Hll region NGC 3603. The effects of field contamination, individual 
reddening, and a possible binary contribution are investigated. The MF slopes 
resulting from the different methods are compared, and lead to a surprisingly 
consistent cluster MF with a slope of F = —0.9 ± 0.15. Analyzing different 
radial annuli around the cluster core, no significant change in the slope of the 
MF is observed. However, mass segregation in the cluster is evidenced by the 
increasing depletion of the high-mass tail of the stellar mass distribution with 
increasing radius. We discuss the indications of mass segregation with respect to 
the changes observed in the binned and cumulative stellar mass functions, and 
argue that the cumulative function as well as the fraction of high- to low-mass 
stars provide better indicators for mass segregation than the MF slope alone. 
Finally, the observed mass function and starburst morphology of NGC 3603 YC is 
discussed in the context of massive local star-forming regions such as the Galactic 
Center Arches cluster, R 136/30 Dor in the LMC, and the Orion Trapezium 
cluster, all providing resolved templates for extragalactic star formation. Despite 
the similarity in the observed MF slopes, dynamical considerations suggest that 
the starburst clusters do not form gravitationally bound systems over a Hubble 
time. Both the environment (gravitational potential of the Milky Way) and 
the concentration of stars in the cluster core determine the dynamical stability 
of a dense star cluster, such that the long-term evolution of a starburst is not 
exclusively determined by the stellar evolution of its members, as frequently 
assumed for globular cluster systems. 

Subject headings: HII regions: individual (NGC 3603) — open clusters and associa- 
tions: individual (NGC 3603, HD 97950) — stars: pre-main-sequence - stars: mass 



1. Introduction 

High-mass star formation is still a puz- 
zle to theorists and observers alike. Star- 
formation theory suggests that high-mass 
stars (M > 20 Mq) cannot form in the 
standard fragmentation and subsequent 
accretion scenario without serious modi- 
fications to the physical processes involved 
(Stahler et al. 2000). Competitive accre- 
tion (Bonnell et al. 2001), merging of pro- 
tostellar clumps or young stars (Bonnell et 
al. 1998), and enhanced accretion in de- 
pendence on the star-forming environment 
(Behrend & Maeder 2001) are all scenar- 
ios suggested for the formation of 0-type 
stars. All these processes share the re- 
quirement of exceptionally high gas and/or 
stellar densities to create high-mass stars. 

The necessity of a high-density envi- 
ronment where altered physical processes 
shape the resultant stellar mass, and pre- 
dominantly create high-mass stars at the 
cost of low- or intermediate mass objects, 
predicts primordial mass segregation in the 
densest star-forming loci, which might be 
evidenced in a flattened initial mass func- 
tion (IMF) with a strong bias to high- 
mass stars in comparison to more mod- 
erate sites of star formation. In contrast 
to this, the IMF is found to be surpris- 
ingly invariant in a diversity of Milky Way 
and Magellanic Cloud star forming regions 
(Massey et al. 1995a, b). 

The densest environments available to 
study the stellar mass spectrum are star- 
burst clusters. In the Milky Way, only 
very few compact, dense starburst clusters 
are known, with the youngest and most 
condensed systems being the Arches clus- 
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ter close to the Galactic Center (GC) and 
the central starburst cluster in NGC 3603. 
These clusters are particularly important 
for our understanding of extragalactic star- 
forming regions, where the starburst mode 
is the predominantly observed mode of star 
formation due to the intrinsic brightness 
of compact, massive young clusters. Star- 
forming regions in other galaxies are, how- 
ever, in most cases too distant to be re- 
solved into individual stars. A detailed 
study of starburst clusters in the Milky 
Way provides a resolved template for ex- 
tragalactic star-forming regions. 

With the aim to observe possible devia- 
tions or to conflrm the universality of the 
initial mass function in starburst environ- 
ments, we derive the stellar mass distribu- 
tion in NGC 3603 YC (Young Cluster). We 
analyze possible changes in the MF slope, 
in particular with respect to the distance 
from the cluster core, where the highest 
densities are found. A possible deviation 
from a standard MF in starburst environ- 
ments would have severe consequences on 
our understanding of star formation in the 
nuclei and tidal interaction zones of dis- 
tant galaxies, such as e.g. the Antennae 
galaxies, and star formation in the early 
universe when the star formation efficiency 
was more intense. This latter aspect refers 
not only to the distant universe, but also 
to the formation epoch of globular clusters 
in the Milky Way. 

As the NGC 3603 YC starburst clus- 
ter is extreme among today's Milky Way 
massive star-forming regions, and with a 
stellar cluster mass of ~ 10^ Mq residing 
in a cloud of total gas mass 4 • IO^Mq 
(Grabelsky et al. 1988) at the low-mass 
end of extragalactic starbursts, several at- 



tempts have been made to understand 
the stellar mass function. In all these 
studies, the resolution in the core was 
(and still is) the limiting factor, causing 
large uncertainties in the resultant MF 
shape. From optical photometry, Mof- 
fat et al. (1994) obtained a MF slope of 
r = -1.4 ± 0.6 for 30 < M < 60 M©, 
i.e., the 77 brightest stars. Hofmann et 
al. (1995) derived F = —1.6 in the mass 
range 15 < M < 50 Mq from only 28 cen- 
tral stars resolved with speckle interfer- 
ometry. While these slopes are consistent 
with a Salpeter MF, only the massive stel- 
lar population is traced. The limitations 
in resolution additionally constrain these 
results mostly to the outer cluster regions 
(except for the speckle data, where statis- 
tics are very poor). The earliest attempt 
towards a complete central MF was per- 
formed by Eisenhauer et al. (1998), who 
derived a MF slope of F = —0.73 for 

1 < M < 85 M© in the innermost 13 
of the cluster center from ESO/ADONIS 
high-resolution adaptive optics photom- 
etry. In this study, crowding effects were 
taken into account via spatially varying in- 
completeness ratios, and the presumed age 
spread of 0.3 — 1 Myr on the PMS was con- 
sidered by calculating a median age of the 
young, intermediate mass population of 0.5 
Myr from individual positions of PMS stars 
in the color-magnitude diagram. Although 
these authors give strong arguments for a 
young main-sequence cluster age of only 1- 

2 Myr, they use a 3 Myr Geneva isochrone 
to derive the MF at the high-mass end. 
Recently, Sagar et al. (2001) confirm this 
slope with F = —0.84, however from ill- 
resolved UBVRI data, such that, again, 
the innermost cluster region could not 
be traced. Despite these uncertainties 



in previous studies, high-resolution V RI 
HST/PC data (FWHM < O'.'l) yield a 
similarly fiat slope ofF = — 0.9±0.1 in 
the mass range 1 < M < IQQ Mq ^initial 
in the inner 20 of the cluster center us- 
ing a composite isochrone for the MS and 
PMS population with a single age of 1 Myr 
(Sung & Bessell 2004). 

We will use the cluster characteristics 
determined in Stolte et al. 2004 (hereafter 
Paper I, where we give a detailed intro- 
duction on NGC3603; see also the intro- 
duction in Sung & Bessel 2004 for a re- 
view on the cluster), in particular a uni- 
form cluster age of 1 Myr and a distance 
of 6 kpc (Paper I), to derive mass func- 
tions and possible radial variations in the 
MF from VLT/ISAAC near-infrared data 
(see Paper I, and Brandl et al. 1999) ob- 
tained under excellent seeing conditions 
(FWHM < 0'.'4). This data set contains 
the deepest currently available photome- 
try of NGC3603YC, allowing one to ob- 
serve cluster members down to masses of 
~ 0.1 Mq, well into the subsolar regime. 

In Sec. 2, we give a brief summary of the 
observations presented in Paper I. The use 
of high-resolution near-infrared data allows 
us to derive the MF slope down to 0.4 Mq 
in the crowding-limited cluster field, while 
at the same time minimizing effects from 
variable extinction. We recall the critical 
cluster properties used in the MF deriva- 
tion in Sec. 3 from Paper I. The mass func- 
tions are derived in Sec. 4, and the resul- 
tant slopes are discussed in Sec. 5. Radial 
variations in the MF slope are analyzed 
in Sec. 6, and the effects of mass segre- 
gation are interpreted in Sec. 7. Because 
of the impact of mass segregation on our 
understanding of the star formation pro- 



cess, we present a detailed comparison be- 
tween different star-forming environments 
in Sec. 8. Finally, we summarise our main 
results briefly in Sec. 9. 

2. Observations 

A detailed analysis of the photometry 
of NGC3603YC is presented in Stolte et 
al. (2004, Paper I), where color-magnitude 
and color-color distributions of the clus- 
ter center were derived and discussed. For 
the sake of completeness, a brief sum- 
mary of the observations contributing to 
the derivation of the mass function is given 
here. Deep VLT/ISAAC observations were 
taken in the Jg, H, and Kg bands. The Jg 
and Ks observations resulted in the best 
resolution due to stable observing condi- 
tions, while the seeing degraded and var- 
ied during the if -band observations. The 
Js and Ks data display the highest spatial 
resolution on the final images, thus yield- 
ing the best photometric performance, and 
will therefore form the basis for the mass 
function derivation. In order to avoid arte- 
facts and nebulous peaks, it was neverthe- 
less required that each object in the final 
photometric catalogue be detected in all 
three filters, including H . 

With the goal to resolve the cluster 
center, two sets of images were created. 
As crowding (as opposed to photometric 
depth) is the most limiting factor in the 
cluster center, the lowest seeing images 
were combined to Js and Kg integrations 
of 9 and 11 min, respectively, covering a 
2x2 field of view centered on the clus- 
ter. The final resolution on the combined 
images is 0. 38 in Jg and 0. 35 in Kg. In ad- 
dition, the total set of 35 images in Jg and 
39 images in Kg yielded combined frames 



of 35 and 39 min exposure time and a field 
of view of 3.'4 x 3.'4, with final resolutions 
of 0. 40 and 0. 37, respectively. This larger 
field is used to estimate the field contami- 
nation in the cluster sample from the edge 
of the field of view, where crowding does 
not affect the photometry. The bright- 
est stars in the field and in particular the 
bright core of NGC3603YC inside a ra- 
dius oi R = 7 suffer from severe satura- 
tion, such that the core itself had to be ex- 
cluded from the analysis. As the most mas- 
sive stars reside inside the core, this implies 
that the upper end of the mass function 
cannot be derived from these data. As no 
NIR data set with comparable resolution 
covering the larger cluster area is available 
at the present time, we cannot easily fill in 
the high-mass end. Analysis of the upper 
end of the mass function in the cluster core 
as determined from high-resolution adap- 
tive optics or HST/WFPC2 data is dis- 
cussed elsewhere (Eisenhauer et al. 1998; 
Sung & Bessell 2004), and the results from 
these studies will be included in the discus- 
sion of the mass function. 

The details of the observations can be 
found in Paper I, Sec. 2, Tabs. 1 and 2. 

3. Cluster properties of NGC3603 YC 

In Sec. 4, various methods will be used 
to derive the mass function in NGC 3603 YC. 
Before the methods suggested from the 
complex physical properties observed in 
the cluster center are described (cf. Paper 
I), the main physical parameters such as 
the distance, the age and the metallicity of 
NGC 3603 YC are recalled from Paper I. 



3.1. Distance, extinction and age 
derived from the PMS transi- 
tion region 

The distance, extinction and age of the 
starburst cluster were determined from 
the ISAAC data via isochrone fitting 
to the pre-main sequence/main sequence 
(PMS/MS) transition region in Paper I. 
Palla & Stahler (1999) pre-main sequence 
isochrones of ages 0.3-3 Myr were used to 
perform the fit. The selection of isochrones 
is discussed in Paper I. As these param- 
eters are crucial for the mass function 
derivation, we summarise the results be- 
low. 

3.1.1. Distance to NGC 3603 

Kinematic distance estimates to NGC 3603 
range from 6 to 10 kpc, while early pho- 
tometric studies yielded ~ 7 ± 0.5 kpc 
(Moffat 1983; Melnick et al. 1989). The 
purely empirical distance of 7 kpc derived 
by Moffat (1983) from optical photome- 
try and spectral classification of massive 
cluster members was essentially used in 
all subsequent studies. Although indepen- 
dent of stellar evolution models, these early 
photometric distance determinations were 
restricted to the main sequence popula- 
tion. 

We obtain a distance modulus of DM = 
13.9 ±0.1 mag from isochrone fitting of 
the PMS/MS transition region. The cor- 
responding distance, rf = 6 ± 0.3 kpc, is 
slightly lower than previously derived val- 
ues, but consistent with the kinematic dis- 
tance estimate of 6.1 ± 0.6 kpc from De 
Pree et al. (1999), and the photometric dis- 
tance of 6.3 ± 0.6 kpc derived by Pandey 
et al. (2000). 



3.1.2. Foreground extinction to NGC 
3603 

The foreground extinction as derived 
from spectra of high-mass stars close to 
the cluster center is Ay = 4.5 ± 0.3 mag 
(Moffat 1983), and Ay = 4.6 ± 0.6 from 
optical photometry (Melnick et al. 1989). 
Melnick et al. 1989 observe a radial ex- 
tinction variation, confirmed by Pandey et 
al. 2000, which suggests an extinction of 
4 mag {Eb^v ~ 1-3, albeit uncertain due 
to large scatter) when extrapolated to the 
cluster center^. 

A foreground extinction of Ay = 4.5 ± 
0.6 mag is derived from PMS isochrone fit- 
ting towards the pre-main sequence pop- 
ulation in Paper I, consistent with ear- 
YC lier studies. The large uncertainty mir- 
rors the low sensitivity of infrared obser- 
vations to granular material along the line 
of sight. This choice also fits the lower 
PMS stars towards fainter magnitudes, be- 
yond the distinct, horizontal transition re- 
gion. A slightly lower extinction of Ay = 4 
mag is derived for main sequence stars, 
pre-dominantly found close to the cluster 
core, slightly lower than extinction esti- 
mates from spectroscopic studies, where 
values around Ay = 4.5 mag were sug- 
gested (Moffat 1983), but consistent with 
the core extinction expected from the ra- 
dial trend observed by Melnick et al. 1989. 
In this estimate the assumption was made 
that the resolved secondary sequence is 
comprised of a physically distinct popula- 



^A radial variation may be an oversimplification, 
especially because IRAS maps indicate a North- 
South orientation of dense material around the 
cluster (Niirnberger et al. 2002). The insensitiv- 
ity of NIR data prohibits a more accurate deter- 
mination of the spatial extinction variation. 



tion (due to binarity or other physical pro- 
cesses), and therefore excluded from the 
isochrone fit. Objects on the secondary 
sequence display redder colors, such that 
the extinction estimate increases when the 
whole population is considered, in agree- 
ment with earlier results. ^ As stars will be 
dereddened individually in the MF deriva- 
tion, which implies shifting stars along the 
reddening path to Ay = mag, the differ- 
ence in extinction between PMS and MS is 
meaningless for the mass calculation. 

3.1.3. Age of the starhurst cluster 

The age of NGC3603YC has evolved 
in the literature during the past 20 years. 
Early studies suggested an age of 2-3 Myr 
from the fact that presumably evolved 
stars with strong winds and Wolf-Rayet 
characteristics reside in the cluster center 
(Moffat 1983). In the past decade, how- 
ever, it became increasingly clear that ex- 
tremely massive stars of type 02 show 
strong Wolf-Rayet characteristics while 
still on the main sequence (see, e.g., Wal- 
born et al. 2002). A main sequence nature 
is evidenced in the spectra of the most 
massive stars in NGC3603YC by strong 
hydrogen emission (Drissen et al. 1995; 
Drissen 1999). Moffat et al. (2002) ar- 
gue that the WR stars in NGC3603YC 
are "probably main sequence stars of ex- 
tremely high luminosity with strong wind- 
produced emission lines" from the cor- 



^Note, however, that increasing color terms for red- 
der PMS stars could not be taken account due to 
the limited number of standard stars. Thus, we 
cannot entirely rule out an instrumental color ef- 
fect between MS and PMS, or uncertainties in the 
extinction law or PMS isochrones, to cause the de- 
rived deviation in Ay. See Paper I for a detailed 
discussion of the uncertainties. 



relation observed between X-ray and ra- 
dio flux, which is uncorrelated for normal 
Galactic (single) WRs. Massey & Hunter 
(1998) compare the spectra of three WRs 
in NGC3603YC to similar stars observed 
in R 136/30 Dor, and conclude that these 
are probably "super O stars". From the 
short main sequence lifetime of these mas- 
sive beasts, the estimated cluster age has 
decreased in the literature to about 1 Myr 
(Drissen 1999; Moffat et al. 2002), which 
is consistent with age estimates indicated 
by PMS isochrones (Brandl et al. 1999; 
Eisenhauer et al. 1998). 

In Paper I, we obtain an age of 1 Myr 
as the best fit to the PMS/MS transition 
region. The significant width of the tran- 
sition region caused the suggestion of an 
age spread between 0.3 and 1 Myr in the 
starburst cluster in these previous studies. 
As discussed in detail in Paper I, there 
are two straight-forward interpretations of 
the width of the transition region. In ad- 
dition to an age spread, the PMS transi- 
tion can also be interpreted as comprised 
of a mixture of single and binary stars, 
with the well-populated upper envelope in- 
dicating binaries with mass ratios close to 
unity. This interpretation is supported by 
a secondary sequence observed close to the 
main sequence in the ISAAC data. From 
X^ distance minimisation, the best fit to 
the PMS population is obtained for a single 
cluster age of 1 Myr, including the assump- 
tion that stars on the secondary sequence 
are binaries. An age distribution can, how- 
ever, not be entirely excluded prior to a 
detailed spectral analysis of objects in this 
regime. For the mass function derivation, 
we will therefore assume a single cluster 
age of 1 Myr, while treating the secondary 



sequence in different ways to minimise a 
potential bias introduced by the interpre- 
tation of the secondary sequence, e.g. as 
binaries. 

As stars on the MS do not evolve 
severely, slight deviations in the age deter- 
mination are not critical for the derivation 
of the main sequence MF (M > 4 Mq ac- 
cording to a 1 Myr Palla & Stabler (1999) 
isochrone). At these young ages, stellar 
evolution is almost negligible except for 
the highest mass objects, which are not 
included here, such that the mass func- 
tion derived from a 1 Myr MS isochrone 
should not deviate significantly from ear- 
lier MFs using a 3 Myr model. The main 
sequence portion of the MFs derived in 
the following is therefore comparable to 
the results of previous studies, while the 
pre-main sequence MF depends strongly 
on the assumed age and the treatment of 
stars in the transition region. 

The isochrones used to determine masses 
of individual stars from the color- magnitude 
plane will be shown along with the PMS 
and MS selection of stars in Sec. 4, Fig. 5 
(see also Paper I). 

3.14. MetalUcity of NGC 3603 

No thorough analysis of NGC 3603 with 
respect to metallicity is available in the lit- 
erature. The spectral analysis of the WR- 
like massive component in the cluster core 
suggests a metallicity close or equal to so- 
lar (Schmutz & Drissen 1999). NGC 3603 
is located at a Galactoccentric distance of 
8 kpc, comparable to the Sun. The ra- 
dial gradient in metallicity observed in the 
Milky Way supports that the metallicity in 
NGC 3603 is comparable to the solar value, 
which was adopted. 



3.2. Uncertainties in the mass func- 
tion derivation 

The color and brightness of stars at 
the youngest evolutionary stages are influ- 
enced by a diversity of processes. Vari- 
able extinction can be caused by circum- 
stellar envelopes or a variation in the den- 
sity of remnant molecular material along 
the line of sight. Near-infrared excess emis- 
sion from circumstellar disks alters the col- 
ors and increases the apparent reddening 
of stars, thereby causing an overestimate 
in the stellar mass of individually dered- 
dened objects with intrinsic NIR excesses. 
Binarity/multiplicity increases the bright- 
ness, again resulting in an overestimate of 
the stellar mass. The large scatter ob- 
served in particular at fainter magnitudes 
in the PMS/MS transition and PMS region 
of the CMD suggests that these evolution- 
ary stages are especially affected, although 
the increasing photometric uncertainty to- 
ward fainter magnitudes prohibits a quan- 
titative analysis of increasing NIR excesses 
in PMS vs. MS stars. These uncertainties 
for masses below 4 Mq have to be kept in 
mind during the derivation and interpreta- 
tion of mass functions. 

While variable extinction can be taken 
into account by individually dereddening 
stars to the theoretical isochrone, NIR ex- 
cess and binarity are difficult to assess. As 
the majority of stars in NGC 3603 YC are 
found close to the fitted isochrone, severely 
reddened sources, either excess cluster 
members or background objects, can be 
excluded from the MF derivation without 
significant statistical loss. Multiplicity can 
usually not be taken into account in MF 
derivations. With the interpretation of 
the secondary sequence in NGC 3603 YC 



as a sequence of binaries with mass ratios 
close to unity, we can, however, estimate 
the contribution of (equal-mass) binaries 
to the MF and correct for this effect. 

4. The Mass Function in the central 
NGC 3603 cluster 

4.1. General remarks on the mass 

function derivation in NGC 3603 YC 

As a consequence of the clear separa- 
tion of a pre-main sequence and a main 
sequence population observed in the cen- 
tral CMD of NGC 3603 YC, and the exis- 
tence of a visible secondary sequence inter- 
preted as binary candidates (see Paper I), 
the mass function derivation has to incor- 
porate a number of physical assumptions. 
While generally one favourite set of stellar 
evolution models is chosen to derive the 
mass function via isochrone fitting, in the 
case of NGC 3603 YC two different sets of 
evolutionary models are required to repre- 
sent the MS and the PMS population at 
a given age. As the 1 Myr pre-main se- 
quence isochrone from the Palla & Stahler 
set of models fits both the transition of 
the PMS to the MS as well as the lower 
PMS distribution best (cf. Paper I), this 
isochrone is chosen for the transformation 
of luminosities into masses on the PMS. 
For the evidence given above and in Pa- 
per I, the MS population is also consistent 
with a single star-forming burst 1 Myr ago. 
Along with the binary interpretation, this 
suggests a single age for the entire clus- 
ter population, without the need to em- 
ploy an age spread. The mass-luminosity 
relation for main sequence stars is deduced 
from a 1 Myr solar metallicity isochrone 
from the Geneva basic grid of stellar evo- 



lution models (Lejeune & Schaerer 2001). 
Effects such as enhanced mass loss are not 
taken into account, as the ISAAC data 
are saturated for stars with masses above 
20 Mq, below which mass loss is not signif- 
icant at these young ages. In addition, the 
case of a 2 and 3 Myr MS isochrone, used in 
prior MF derivations and allowing for some 
uncertainty in the derived age, will also be 
tested to deduce potential effects of stellar 
evolution on the resulting MF slope. 

The identification of a secondary se- 
quence poses an additional challenge to 
a realistic MF derivation. As no prior 
data set of the starburst cluster resolved 
this feature clearly, the MF will be shown 
with and without binary correction to al- 
low comparison with prior MF derivations. 

As a consequence of the variable extinc- 
tion observed in the NGC 3603 region, in- 
dividual dereddening was applied to deter- 
mine stellar masses. As the NIR data alone 
do not allow us to derive the extinction law 
toward NGC 3603, the extinction law was 
adopted from Rieke & Lebofsky (1985). 
Each star has been shifted along the red- 
dening vector in the Jg, Jg — Kg plane onto 
the PMS or MS isochrone. Although this 
procedure is more accurate in the sense 
that individual reddening is taken into ac- 
count, it implicitly assumes that all scatter 
observed in the CMD is caused by redden- 
ing variation. This is clearly not the case. 
Besides the photometric scatter inherent in 
the data, other processes such as binarity, 
stellar rotation, and intrinsic IR excesses 
contribute to the scatter. Most of the ad- 
ditional processes influencing a very young 
stellar population are difficult to quantify, 
and cannot be taken into account in the 
following analysis. While photometric un- 



certainties are expected to introduce ran- 
dom scatter, such that their combined ef- 
fect should be neghgible, individual dered- 
dening allows one to estimate the maxi- 
mum influence of extinction variations on 
the MF slope. The intense crowding in the 
cluster center decreases the detectability of 
lower-mass stars and can cause a system- 
atic uncertainty in the MF slope, which 
would mimic a bias to high-mass stars. For 
each cluster area/annulus the average in- 
completeness is therefore estimated as a 
function of the radial distance from the 
cluster center and hence stellar density. 

While a detailed analysis of the uncer- 
tainty in the MF slope imposed by the 
choice of stellar evolutionary models is be- 
yond the scope of this paper, a thorough 
discussion on this subject can be found in 
Andersen (2004). 

The procedures used to derive the stel- 
lar mass distribution in NGC 3603 YC are 
outlined in detail in the next sections. 

4.2. Incompleteness correction 

Artificial star tests were performed on 
each science frame to determine the recov- 
ery fraction per magnitude. Between 200 
and 250 stars were inserted in each artifi- 
cial frame to avoid enhancing crowding ef- 
fects. 45 artificial frames were created from 
the smaller 2x2 "cluster" images and 30 
artificial frames from the larger 3.4 x 3.'4 
"field" images. The inserted number of 
stars corresponds to less than 1% of the 
total number of stars detected on each im- 
age. Magnitudes and positions were as- 
signed randomly using the addstar task in 
IRAF. As the cluster itself covered only a 
small area on the high-resolution frames, 
additional simulations were performed on 



the innermost 60 x 60 , with 100 to 150 
artificial stars or < 2% of the stellar pop- 
ulation found in the cluster center. Source 
detection and PSF fitting was conducted 
with the same parameters as for the stellar 
photometry. The procedure was repeated 
to yield the recovery fractions for ~ 6000 
stars in the field frames, and ~ 11000 stars 
in total on the cluster frames, with 6000 
stars distributed randomly over the entire 
image, and 5000 stars confined to the crit- 
ical central region. The resultant recovery 
fractions were applied in two steps during 
the analysis: 1) during field star correction 
and 2) during incompleteness correction of 
the mass function. 

4.3. Field star correction 

The field star fraction was derived from 
an area to the North and East of the clus- 
ter center. This region is particularly well 
suited to estimate field stars in the low- 
extinction cluster center, as it is neither 
infiuenced by strong Hll region emission 
nor shows evidence of variable or enhanced 
extinction. Field stars were statistically 
subtracted in each 0.5 x 0.5 mag color- 
magnitude bin of the CMD. Prior to sta- 
tistical field star subtraction, the number 
of field stars had to be corrected for the 
varying stellar density (crowding) with ra- 
dial distance from the cluster center. For 
this purpose, the field and cluster recovery 
fractions, f field and f duster, were fitted as 
a function of the magnitude as a third or- 
der polynomial as displayed in Fig. 3. The 
number of field stars is corrected for field 
incompleteness as nfieid,corr = nfieid/ f field, 
i.e. to a field completeness of 100%. 
This number was then reduced accord- 
ing to the expected recovery fraction of 
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stars in the radial annulus studied, yield- 
ing a combined correction of Ufieidjinai = 

J cluster ' '^field,corr J cluster / J field ' ''' field) 

with ffieidjinai = f duster / f field corresponds 
to the final density/crowding correction 
factor applied to adjust the field star num- 
ber counts to the cluster center. In addi- 
tion, the measured number of field stars 
on the larger field area was scaled down to 
the cluster or annulus area considered. 

The incompleteness simulations yielded 
independent correction factors for each fil- 
ter. In principle, the combined probability 
to find a field star in each color-magnitude 
bin in the Js, Js — Kg plane is then given 
by the product of the individual probabil- 
ities, fjK = fj ■ fx- In the color and mag- 
nitude range covered by the cluster cen- 
ter stars, however, the incompleteness is 
entirely determined by the J-band limit. 
Thus, only J-band incompleteness is taken 
into account in the mass function deriva- 
tion. The correction factor fj for each 
magnitude is obtained from the polynomial 
fit. The measured number of field stars 
is corrected by the factor derived for each 
annulus studied, and the adjusted number 
of field stars is then statistically removed 
from the cluster CMD. 

As the extent of the cluster is de- 
bated, with estimates ranging up to a ra- 
dius of 2.5 (Niirnberger & Petr-Gotzens 
2002), the area covered by the ISAAC data 
may overestimate the field star contami- 
nation. As the low-mass halo surround- 
ing NGC3603YC is mainly composed of 
faint, low-mass stars, the bright central 
cluster population should not be severely 
affected^. 



Fig. 4 shows the effects of field subtrac- 
tion on the central cluster population, and 
also indicates the 50% completeness limit 
for this radial selection, down to which 
mass function slopes will be derived. 

In the following sections, the different 
approaches used to derive the mass func- 
tion in NGC3603YC will be described in 
detail. Results for all MF derivations are 
summarised in Tab. 2. The combined 
PMS/MS mass functions of this selection 
are fitted in Fig. 6. The apparent MS pop- 
ulation fainter than Jg = 15.5 mag, cor- 
responding to the 1 Myr MS turn-on at 
~ 4:Mq, is likely comprised of MS stars 
in the NGC3603 region. This population 
has to be far older than the cluster age, 
and cannot originate in the cluster. The 
lower MS beyond the cluster MS turn-on 
is therefore excluded from the fit. Thus, 
the combined low-mass MF below the MS 
turn-on at M = 4Mq is solely comprised 
of PMS stars. 

4.4. Simple star counts 

In the simplest model, the luminosity 
of a star refiects its evolutionary stage and 
mass. This ignores all processes chang- 
ing the luminosity and thus magnitude of 
a star individually and independently of 
its evolutionary stage. The isochrone can 
then be applied directly to each given stel- 
lar magnitude to transform the brightness 
distribution into a mass function. As this 



■^Ideally, a separated field should be used to es- 
timate the background population. One ofF-field 



was observed at a distance of ~ 30 from the clus- 
ter center, but unfortunately turned out to contain 
a different stellar population dominated by older 
giant stars. The strong spatial variations in the 
stellar population along the Carina arm compli- 
cate the determination of the field star contribu- 
tion in any particular region. 
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is the most frequently used method to de- 
rive mass functions, we include it here as 
the simplest approach. Field data are not 
always available, and field contribution can 
indeed be negligible in the dense cluster 
regions, such that field star contamination 
has rarely been taken into account in MF 
derivations of the NGC3603YC (except 
for the study of the cluster center by Eisen- 
hauer et al. 1998). In order to show the 
effects of field contamination alone, MFs 
were calculated from simple star counts 
both with and without field subtraction. 



4.5. Individual dereddening 

The stellar population in NGC3603 
shows individually varying extinction, caused 
either by variations in the distribution 
of granular material in the foreground, 
or by individual reddening from remnant 
circumstellar envelopes. The maximum 
possible effect of individual reddening on 
the MF slope can be probed by assum- 
ing that the photometric scatter is caused 
by varying extinction alone. Individual 
reddening is then taken into account via 
dereddening each star along the redden- 
ing vector (Rieke & Lebofsky 1985) in the 
Js,Js — Kg plane until an intersection with 
the corresponding isochrone is reached, 
and the mass of the star is determined 
from the intersection point. The extinc- 
tion applied to the isochrone itself is not 
relevant in this case, as only the length 
of the dereddening path (i.e., the relative 
shift in the CMD) depends on the location 
of the isochrone, but not the desired in- 
tersection point. For the MF derivations 
below, isochrones matching the average 
foreground extinction of the MS and PMS 



populations were used, and stars blueward 
and redward of the isochrone have been 
shifted accordingly. The blue and red lim- 
its in the distance from the isochrone was 
determined by eye from the stellar distri- 
bution in the Js-,Js — Kg plane. As the 
main sequence is not affected by field star 
contamination, all stars around the MS are 
included in the MF, with a color selection 
of -0.25 mag blueward and -1-0.35 mag red- 
ward of the MS isochrone. The bulk of the 
PMS population, where reddening is more 
severe, can be included selecting colors - 
0.25 mag blueward and -1-0.5 mag redward 
of the PMS isochrone. This rejects the 
reddest objects {Ay > 7.5 mag), which 
might be affected by infrared excesses and 
background contamination. The selection 
of stars contributing to the mass function 
is shown in Fig. 5. 

Note that most of the stars are found in 
a very narrow range around the PMS, in an 
envelope of only -0.1 to -1-0.15 mag. This 
narrower selection does not alter the slope 
of the PMS MF alone, but in comparison 
with the MS MF too few stars are found 
in the transition region, such that MS and 
PMS do not agree at the transition mass. 
Therefore, the wider PMS selection was 
used, resulting in very good agreement of 
the MS and PMS MF bins at the transi- 
tion mass. 

In all following MF derivations, field 
subtraction and individual reddening are 
taken into account. 

4.6. Treatment of the PMS/MS 
transtion region 

Stars in the transition region between 
PMS and MS were allocated PMS transi- 
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tion masses, which are in the narrow range 
between 2.5 and 4:Mq, depending on their 
distance to the isochrone. No individual 
dereddening was performed in this case due 
to the large scatter in the transition re- 
gion. There is a slight ambiguity in the 
Palla & Stabler isochrone in the sense that 
stars with 3.5 to 4Mq cannot be distin- 
guished, but this is not relevant for the re- 
sultant MF as these stars all fall into the 
same mass bin. Only masses above the hy- 
drogen burning turn-on of 4 Mq enter the 
MS MF. Although this facilitates the anal- 
ysis as a clear cut is introduced between 
MS and PMS contributions, the transition 
region is clearly physically more complex 
than envisioned in this simple treatment. 
In particular, stars contributing to the MS 
may be scattered away from the narrow 
range beyond the PMS transition point in 
the CMD, rendering the respective MS and 
PMS contributions close to the 4 Mq PMS 
limit very uncertain. The transition mass 
will be marked in all derived MFs, such 
that the combination of MS and PMS con- 
tributions is clear, and in particular the 
last PMS and first MS bin should be ob- 
served with care. 

4.7. Treatment of the secondary se- 
quence 

Evidence for a binary nature of a sec- 
ondary sequence observed in the ISAAC 
CMD of NGC3603YC is given in Paper 
I. However, this interpretation is as yet 
unproven by spectroscopy. Three differ- 
ent approaches are therefore chosen in the 
treatment of the secondary sequence. 1) 
The secondary sequence is completely re- 
jected from the mass function derivation. 
As the relative fraction of stars on the sec- 



ondary sequence is not observed to vary 
with stellar brightness within the statis- 
tical uncertainties, removing stars on the 
secondary sequence corresponds to statisti- 
cally subtracting ~ 30% of the population 
above 2.5 Mq, below which the vertical 
turn of the PMS does not allow one to dis- 
tinguish a secondary sequence. While this 
may result in too low star counts in high- 
mass bins with respect to the low-mass 
MF, the statistical properties and thus the 
slope in the high-mass MF should not be 
altered. 2) The distinct interpretation is 
ignored and the entire population is used 
for the MF derivation, in agreement with 
earlier studies. This MF is also directly 
comparable with MF derivations that do 
not have the photometric accuracy or reso- 
lution to attempt a binary correction, such 
as the early study by Salpeter (1955). 3) 
The binary interpretation is used to cor- 
rect for the effect of binaries with mass 
ratios close to unity on the resultant MF 
slope. 



4-7.1. Binary rejection 

The selection of stars on the secondary 
sequence was shown and carefully dis- 
cussed in Paper I, see esp. Fig. 7. The band 
of stars above the transition region and 
redward of the main sequence was rejected 
in the first estimate of the MF, assum- 
ing that some currently unknown physical 
process causes the offset in the CMD, and 
thus masses may not be accurately deter- 
mined. If stars on this sequence are not 
subject to a mass bias from a different for- 
mation process, the resultant "single" mass 
function should represent the underlying 
statistically-unbiased stellar population of 
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the star cluster. The vertical distribution 
of lower PMS stars does not allow such a 
selection. Even if the same photometric 
offset is included in the case of PMS stars, 
as is to be expected from the secondary 
sequence extending all the way into the 
transition region, such an offset is invis- 
ible in the vertical PMS feature. As a 
consequence, the corresponding PMS stars 
cannot be rejected, and only the MS and 
the transition region (down to 2.5 Mq) are 
affected by binary rejection. 

4-7.2. Binary correction - MS and tran- 
sition region 

As a first step toward a binary corrected 
MF, the visible binary candidate sequence 
above the main sequence and the transi- 
tion population was adjusted. The J^-band 
magnitude of the selected stars was in- 
creased by 0.75 mag as expected for equal- 
mass binaries, and an additional compan- 
ion star with equal color and magnitude 
was included in the number counts. This 
procedure makes the simplifying assump- 
tion that all binaries contributing to the 
brighter secondary sequence are close to 
equal-mass systems as suggested from the 
observed magnitude offset, in accordance 
with the observational restriction that for 
arbitrary mass ratios no prediction can be 
made. Thus, the MFs derived below are 
still system MFs uncorrected for compan- 
ions with masses significantly lower than 
the primary mass. The adjustment is ap- 
plied after field subtraction, before individ- 
ual dereddening and MF calculation. This 
way, the bias from individual dereddening 
of equal-mass binaries due to the larger 
brightness and thus higher mass estimate 
is avoided. The resultant pairs of "single" 



stars are treated in exactly the same way 
as the rest of the population. 

4. 7.3. Binary correction - statistical PMS 
correction 

While a secondary sequence is observed 
above the MS and the PMS transition re- 
gion, the problem of the merging of bi- 
naries into the lower PMS population re- 
mains. The fraction of stars observed on 
the secondary sequence is approximately 
30% of the total number of stars both par- 
allel to the MS and in the PMS/MS tran- 
sition region (J < 15.5 mag). When inter- 
vals of 1 mag width are studied, no trend 
is observed (Tab. 1) toward fainter stars. 
Thus, to correct for the equal-mass binary 
contribution on the PMS a constant bi- 
nary fraction of 30% has been assumed. If 
no binary formation mechanism enhances 
the number of high-mass binaries or de- 
creases the number of low-mass binaries 
with mass ratios close to unity, a signifi- 
cant fraction of binaries should be hidden 
in the PMS distribution. In order to cor- 
rect for the bias imposed on the derived 
MF, 30% of the PMS stars were selected 
randomly, shifted down by +0.75 mag and 
an additional companion was added, as in 
the case of the MS binary sequence. If 
the equal-mass binary fraction is the same 
for each mass bin, as suggested from the 
MS and transition region, this procedure 
should statistically correct for the overes- 
timate of individual masses due to brighter 
magnitudes of potential binaries. 

5. Discussion of the resultant mass 
functions 

Mass functions were created from the 
above procedures with a histogram bin 
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width of log(M/M0) = 0.2, starting at 
0.1 Mq. a linear least-squares fit was ap- 
plied to the incompleteness corrected MF 
from the 50% completeness limit all the 
way to the highest mass bin. At the high- 
mass end of the MS MF, saturation losses 
thin out the CMD. As a consequence of the 
combination of different exposures taken 
under varying conditions, the saturation 
limit is not a sharp point, but extends over 
a range of ~ 1 mag in J^. The approximate 
fraction of stars lost due to saturation in 
the upper mass bin, M < 20 Mq, is ~ 50%. 
The high-mass end of the MF has been cor- 
rected accordingly, and is included in the 
MF fit. The MF in NGC3603YC is con- 
sistent with a single-slope power law over 
the entire mass range probed by our data, 
0.4 < M < 20 M©. Note that the apparent 
turn-over below OAMq is not included in 
the fit, and we are not able to draw con- 
clusions on the possible turn-over mass in 
NGC3603YC, because the MF is heavily 
depleted by incompleteness, and the cor- 
rection increasingly uncertain. 

All resultant MF slopes are summarised 
in Tab. 2 and sample MFs are shown in 
Fig. 6. The derived individual slopes of the 
MS alone (see Fig. 6 and Tab. 2) are meant 
to illustrate effects of the different correc- 
tions introduced, and should not be taken 
at face value due to the low number of bins 
included in these fits. The PMS isochrone 
covers a large range in color over a narrow 
mass range in the transition region where 
PMS stars equilibrate to main sequence hy- 
drogen burning. The colors in the transi- 
tion region range from J — K = 0.7 mag to 
J — K = 1.3 mag with masses between 
2.5 and 4Mq. Only four model points 
are available for the transition region, such 



that the mass transformation had to rely 
on interpolation between theoretically pre- 
dicted color, magnitude and thus mass val- 
ues. This is likely the cause for the en- 
hancement of stars seen in each MF in the 
mass bin centered on 2 Mq, and the subse- 
quent depletion in the number of stars in 
the following mass bin centered on 3.2 Mq. 
The combined MS/PMS mass function is 
less affected by the two subsequent bins at 
PMS/MS transition masses, such that the 
continuous mass function yields the most 
reliable fit. 

After the effort undertaken to correct 
for the effects discussed above, the de- 
rived mass function slopes show remark- 
ably small scatter. A few general tenden- 
cies are observed. First of all, the slope 
of the PMS MF is very stable with values 
ranging from P = —0.61 to —0.85, with 
a scatter smaller than the formal least- 
squares fitting error of ~ 0.25. The MS MF 
is more severely affected by binary correc- 
tion than the PMS MF. This is expected 
as the upper mass bins are emptied with- 
out substitution from higher mass bins due 
to the saturation truncation, while binary 
intermediate- and low-mass stars are re- 
plenished from brighter bins. As the up- 
per end of the MF is heavily affected by 
saturation, the loss in stars and resultant 
steepening of the MF is a consequence of 
saturating the binary sequence. 

The slope of the combined mass func- 
tion is not considerably affected by the MS 
variations, but shows remarkable consis- 
tency with values ranging from -0.84 to - 
0.99 with a formal fitting error of ±0.1. As 
expected for the core region with R < 1 
pc, field contamination has only a minor 
effect on the MF slope (Tab. 2, cf. Fig. 2). 
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The use of an older main sequence of 2 
or 3 Myr does not change the combined 
MF fit significantly (see Tab. 2), as ex- 
pected^. The most reliable derivation of 
the mass function in the central cluster 
area of NGC3603 YC including all correc- 
tions discussed in the previous sections is 
displayed in the bottom panel of Fig. 6. 
The slope of the mass function including 
field subtraction, individual dereddening 
and binary correction is F = —0.91 ± 0.08 
over the mass range 0.4 < M < 20 Mq. 

A slope of F = -0.77 down to M > 
1 Mq was derived by Eisenhauer et al. (1998) 
from high-resolution ADONIS data in the 
innermost 13 of the cluster center, in 
reasonable agreement with our estimate. 
These authors did not account for binaries, 
and interpreted the observed scatter in the 
PMS/MS transition as an age spread in the 
range 0.3 to 1 Myr. Accordingly, a mean 
age of 0.5 Myr was used to derive masses 
of PMS stars. As shown in Paper I, Fig. 4, 
this interpretation is consistent with the 
data as long as the photometric accuracy 
does not allow one to resolve the equal- 
mass binary offset. A younger age and 
thus evolutionary state results in assigning 
lower masses to the stars. Comparably, bi- 
narity correction results in assigning stars 
with a given brightness a (well-defined) 
lower mass. As the 0.3 Myr isochrone co- 
incides well with the binary sequence, the 
masses on the PMS derived for the tran- 
sition stars are comparable. On the other 



"^Note that for the PMS only the best-fitting model 
was used, while other models and ages unable to 
reproduce the stellar population in the CMD are 
not considered (see Andersen 2004 for a discus- 
sion on MF slope variations as a function of PMS 
models). 



hand, the effect of possible companions on 
the MF is not taken into account in the 
case of an age spread. Statistically adding 
companions at a constant fraction in all 
magnitude bins should, however, not alter 
the MF slope. Thus a comparable slope in 
the MF using a younger isochrone for the 
transformation of magnitudes into masses 
as opposed to binary correction can be un- 
derstood from the similarity of the 0.3 Myr 
and the binary isochrone. 

A severe difference between both data 
sets remains, limiting the above compar- 
ison. Mass segregation in the cluster, in 
the sense that massive stars may be found 
predominantly in the cluster center as ob- 
served in the case of the Arches cluster 
(Stolte 2003), would cause the MF slope 
to steepen with increasing radius. As the 
ISAAC data are limited to R > 7 due to 
saturation, while the ADONIS data only 
extend out to i? < 13 , a steeper MF 
would be expected from ISAAC than from 
the central ADONIS data in the presence 
of mass segregation. The comparison is 
complicated by the different procedures 
used. When the MF is calculated with 
mere field subtraction, without deredden- 
ing and binary correction for a 1 Myr 
isochrone, the slope becomes F = —0.86, 
slightly, but not significantly steeper than 
the ADONIS slope. No significant change 
can be found between the ADONIS core 
MF and the ISAAC MF outside the core. 
This agreement is even more surprising as 
the core is known to be mass segregated 
within R < 6" from HST/WFPC2 data 
(Sung & Bessell 2004). Sung & Bessel de- 
rive a MF slope on the HST/PC area of 
40" X 40" {R < 20") of F = -0.9 ± 0.1 in 
the initial mass range 1 < M < 100 Mq. 
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This slope agrees surprisingly well with our 
slope using present-day masses. Stellar 
evolution at these young ages only becomes 
significant for stars in excess of ~ 50 Mq, 
such that the mass range of 0.1 < M < 
20 Mq covered by our ISAAC data is not 
influenced by post-main sequence evolu- 
tion. The correspondence in both MFs 
might be a cause of the consistent age es- 
timate of a 1 Myr single-age starburst in 
both derivations, but the consistency is 
still surprising in view of the different pro- 
cedures and evolutionary models applied. 
In summary, the integrated MF slope in 
NGC3603YC is with -0.9 slightly flat- 
tened with respect to a Salpeter slope, 
r = —1.35, suggesting the efficient forma- 
tion of high-mass stars in the cluster cen- 
ter, but not a strong deviation from a nor- 
mal IMF. 

6. Radial variation in the mass 
function of NGC 3603 YC 

Radial mass functions were derived in- 
cluding field subtraction, individual dered- 
dening and correction for binaries on the 
MS as well as on the PMS as described 
above. Three annuli covering equal areas 
were selected, each covering 1/3 of the area 
enclosed within 7 < i? < 33 . The small- 
est accessible radius of 7 is determined by 
saturation in the cluster core. The radial 
selections correspond to 7 < i? < 20 , 
20" < R < 27", and 27" < R < 33" (cf. 
Fig. 1). The mass functions for these an- 
nuli are shown in Fig. 7. 

The slopes of the MFs displayed in 
Fig. 7 in the innermost two annuli agree 
roughly within the uncertainties, and only 
a marginal increase from —0.9 ± 0.15 (7 < 
R < 20") to -1.1 ± 0.15 (20 < i? < 27") is 



observed. Sung & Bessel (2004) observe an 
increase from F = — 0.5±0.1 in the cluster 
core {R < 6") to -0.8±0.2 for 6 < i? < 12" 
and -1.2 ± 0.2 for 12 < i? < 20", sug- 
gesting that mass segregation is most pro- 
nounced in the core and the MF steepens 
to normal values rapidly beyond the core 
radius, in good agreement with our values 
for r taking into account the different ra- 
dial selections applied. In the outermost 
annulus, 27 < i? < 33 , the MS is almost 
entirely depleted, and a reasonable power 
law ffi can only be applied to the PMS. 
The PMS ffi yields F = -0.8 ± 0.2, such 
that a further systematic increase in the 
slope at larger radii is not supported by 
our data. 

Despite the comparable slopes in all 
annuli, the truncation and depletion of 
high-mass stars yields strong evidence for 
mass segregation in NGC 3603 YC. This 
indicates that the slope alone is not ad- 
equate to describe mass segregation pro- 
cesses. While the low-mass (PMS) dis- 
tribution decreases in number as the den- 
sity decreases with larger radial distance, 
the shape in each annulus shows a close 
resemblance. The high-mass (main se- 
quence) end on the other hand becomes 
more and more depleted with increasing 
radius. From the incompleteness corrected 
number counts, we can estimate whether 
the MS depletion is consistent with the ex- 
pected decrease in stellar density as ob- 
served for the PMS. While the cluster cen- 
ter (7" < R < 20") harbors ~ 250 PMS 
stars, the PMS population in the second 
annulus decreases to ~ 87 or 35% of the 
central population. In the third annulus, 
72 PMS stars remain, indicating a slow de- 
crease in the density of low-/intermediate- 
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mass stars. The MS population harbors 78 
high-mass stars in the innermost annulus. 
With the same decrease in stellar density, 
27 MS stars are expected in the second an- 
nulus, while only 17 are observed, and 23 
should still be present in the third annulus, 
where only 4 (!) MS stars are found. We 
can conclude that the depletion on the MS 
is much stronger than expected from the 
decrease in stellar density with increasing 
distance from the cluster center, providing 
strong evidence for segregation of massive 
stars inside the cluster core. 

In stark contrast to the strong outwards 
variation in the MS population, the shape 
of the PMS distribution of stars remains 
similar throughout all annuli. The trun- 
cation at very low masses in the first and 
second annulus is due to crowding losses 
of faint stars, as indicated by the incom- 
pleteness correction. As the lowest mass 
bins are incomplete by more than 90%, a 
final conclusion on a possible change in the 
low-mass population can only be drawn for 
stars more massive than IMq, where the 
completeness fraction is more than 50% in 
all annuli. No significant change is ob- 
served in the 3 mass bins above this thresh- 
old. The 50% limit progresses down to 
OAMq in the two outer annuli, where no 
strong variation in the MF is observed be- 
tween both these annuli. Even when the 
analysis is extended to radii of 65 , well 
outside the distinct cluster center, the field 
subtracted MF displays the same slope in 
the PMS regime (see Tab. 3). The rel- 
ative fraction of stars detected down to 
the lowest mass included in the Palla & 
Stabler models, 0.1 M0, increases notice- 
ably in the outermost annulus. Taking into 
account the incompleteness correction in 



all annuli, a truncation in the mass func- 
tion of NGC 3603 YC at the low-mass end 
can safely be ruled out down to 0.4 Mq. 

7. Mass segregation in NGC 3603 YC 

As the slope of the combined MS+PMS 
mass function does not vary significantly 
with radial distance from the cluster cen- 
ter, this measure cannot be used to draw 
conclusions on mass segregation in NGC 3603 YC. 
From the mass functions, there is strong 
evidence that the PMS population alone, 
covering a mass range of 0.4 < M < 4 Mq 
from the 50% completeness limit to the 
PMS/MS transition, is not segregated. 
There is no indication of a flattening or 
depletion of the MF at the low-mass end 
towards the cluster center. Nevertheless, 
the truncation of high-mass MS stars with 
M > 4 Mq beyond a radius of 27 is a 
clear indication of mass segregation in the 
massive population. Such a behaviour can 
be expected from primordial mass segre- 
gation assuming that high-mass stars form 
predominantly in the densest cluster re- 
gion, and may thus have resided close to 
the cluster center throughout their life- 
time, while low-mass stars can form any- 
where inside the parental cloud. How- 
ever, dynamical segregation will also trans- 
port high-mass stars inwards and low-mass 
stars outwards during interactions. 

Although in a cluster as young as ~ 1 
Myr a location of massive stars in the 
cluster center is very suggestive of primor- 
dial segregation, Bonnell & Davies (1998) 
derive a similar stellar distribution from 
N-body simulations of dynamical clus- 
ter evolution. For their richest cluster 
with 1500 stars statistically drawn from 
a combined Salpeter (1955) and Kroupa et 
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al. (1990) MF truncated at miow = 0.1 Mq, 
and a standard number density profile, 
n oc r~^, 90% of tlie massive stars mi- 
grate into tlie innermost 0.5 lialf-mass 
radii witliin 10 crossing times, no matter 
wlietlier tlie massive stars were originally 
randomly distributed in the cluster po- 
tential or concentrated near the half-mass 
radius. Independently of the exact segre- 
gation timescale, these models thus also 
yield strong evidence that initial condi- 
tions are lost rapidly. Although we cannot 
derive the core radius and in particular the 
crucial contribution of the massive compo- 
nent in the core to the cluster potential 
from the ISAAC data, a core radius of 
~ 0.2 pc is estimated from HST/WFPC2 
observations (Grebel et al., private com- 
munication). Furthermore, a velocity dis- 
persion of a few km/s is typically observed 
in dense clusters. This core radius leads to 
a crossing time of only 10^ yr in the case 
of a velocity dispersion of 2 km/s. Ac- 
cording to Bonnell & Davies, heavy mass 
segregation should have created a high- 
mass dominated cluster core after lOt cross 
or in our example only 1 Myr. This order 
of magnitude estimate suggests that the 
segregation timescale expected from clus- 
ter evolution models is comparable to the 
cluster age, and thus that the cluster was 
already shaped by heavy dynamical segre- 
gation. This early loss of initial conditions 
prohibits conclusions on the formation lo- 
cus of high-mass stars. We can there- 
fore not distinguish primordial from dy- 
namical segregation in NGC3603 YC from 
the present-day distribution of high-mass 
stars. 

To derive more quantitative conclusions 
on the mass segregation in NGC3603YC, 



two other means are investigated in the 
following sections. First, cumulative func- 
tions were created from the mass distribu- 
tion to avoid the binning dependence in- 
herent to the standard MF derivation. Sec- 
ondly, the number ratio of high-mass to 
low-mass stars is analyzed as a function of 
the radius. 

7.1. Cumulative Functions 

Cumulative functions (CFs) are cre- 
ated by consecutively adding stars with 
decreasing mass, starting with one star 
at the highest mass observed. The field- 
subtracted, color-selected distribution of 
stars used to derive the mass functions 
enters the cumulative functions. Masses 
are assigned to individual stars during in- 
dividual dereddening in the same way as 
for the MF derivation. The only differ- 
ence to the mass distributions in the MFs 
is that the artificially added companion 
stars for binary correction are not counted 
as stars in the CF. These "stars" have 
not been observed as individual sources, 
and the correction is - in particular on 
the PMS - purely statistical. As the only 
indication for binaries is given by the bi- 
nary candidate sequence observed in the 
higher mass population, only equal-mass 
binaries can be taken into account. Sys- 
tems with smaller mass ratios cannot be 
distinguished from single stars. If the bi- 
nary stars are not biased to a certain mass 
range, the mass distribution in binary stars 
is the same as in single objects. Adding 
companion stars thus adds to the total 
mass, but not to the statistical properties 
of the distribution. The brightness correc- 
tion of 0.75 mag for stars on the secondary 
sequence, however, is crucial, as ~ 30% of 
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the individual masses are overestimated 
if the brightness is indeed enhanced in 
these stars (by binarity or other physical 
mechanims). Thus, the CFs are created 
from all stars selected as cluster members 
in the CMD after field-subtraction, with 
masses corrected for the derived binary 
fraction of 30%. 

The resultant cumulative functions of 
the central cluster population as well as the 
radial variations in the CFs are shown in 
Fig. 8. Theoretical CFs corresponding to 
a single exponent power-law MF are un- 
derlaid for exponents of —0.3 > F > —1.3. 
For the main cluster area, 7 < i? < 33 , 
the main sequence stars in the mass range 
4 < M < 12 M© follow closely the F = 
—0.7 line. The irregular shape at higher 
masses is partially due to random sampling 
of the upper MF as well as saturation losses 
for M > 15M0. Between 3 and AMq, the 
irregularity arises in the ambiguity of stel- 
lar locations in the CMD and thus individ- 
ual masses in the PMS/MS transition re- 
gion. The pre-main sequence distribution 
favours a flatter value of F ~ — 0.5. 

In the cluster center {R < 20 ), the 
CF tangents the F = —0.3 curve for 
M < 2Mq, and F tends to -0.5 at higher 
masses. At larger radii, the CF falls to- 
wards steeper values of F as expected in 
the case of mass segregation. On the 
PMS, however, this effect is weak, and the 
general shape of the CF remains remark- 
ably constant out to radii where almost 
no main sequence stars remain. This is 
particularly pronounced in Fig. 9, where 
the PMS region below 3 Mq can hardly be 
distinguished in all annuli except for the 
slightly flatter CF in the cluster center. 
This supports the flnding of Niirnberger & 



Petr-Gotzens (2002) of a cluster extent of 
-R ~ 150 beyond the edges of our field of 
view. This pre-main sequence population 
can be interpreted as an extended halo of 
low-mass stars comparable to the halo ob- 
served around the Trapezium in the Orion 
Nebula Cluster (Hillenbrand 1997). 

The situation at the high-mass end is 
very different. For stars with masses M > 
2Mq, the CF in the outer annuli indi- 
cates increasingly steeper values of F ~ 
—0.7 to —0.9. At larger radii, the high- 
mass end becomes increasingly depleted, 
evidenced in a steep decline of the CF to- 
ward higher masses. The most striking fea- 
ture in the main sequence population is the 
sharp truncation of the CF at ~ 2.5 Mq for 
radii R > 27 . This cut-off corresponds 
to the absence of main sequence stars in 
the mass functions (cf. Fig. 7), and indi- 
cates that massive stars are highly segre- 
gated toward the cluster center. Beyond 
33 , the irregular CF shape and increase 
in the scaled frequency of high-mass stars, 
obvious in Fig. 9, marks the onset of field 
contamination, as at large radii statistics 
of cluster stars become poorer and realistic 
field subtraction more difficult than in the 
dense cluster center where cluster members 
dominate. 

From these cumulative distributions, we 
see that binning effects influence the slope 
of the mass functions signiflcantly, and the 
assumption of a single exponent power-law 
distribution is over-simplifled. While the 
MF concept is useful as a tool to compare 
results from different data sets and regions, 
yielding a single characteristic value, the 
assumptions entering such analyses have to 
be carefully reviewed before conclusions on 
"universality" or "deviation" from a cer- 
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tain mass distribution are drawn. 

7.2. Fraction of high- to low-mass 
stars 

The fraction of high-mass main se- 
quence stars to low-mass pre-main se- 
quence stars, f (high/low) = n(> 4M0)/n(< 
4:Mq), is derived as a function of distance 
from the cluster center. Individual mass 
derivations for stars in each separate an- 
nulus are used including field-star and in- 
completeness corrections for each radial 
selection. In addition, the cluster popula- 
tion within 7 < i? < 33 , yielding better 
statistics than the outer annuli, is split 
into 5 bins in steps of 2. 5 to allow a more 
complete radial coverage. The MS/PMS 
transition mass, M = 4Mq, as the natural 
mass separation in the stellar population 
of NGC3603YC is chosen to distinguish 
high- and low mass stars. Consequently, 
the fraction of high- to low-mass stars also 
refiects the fraction of MS to PMS stars. 

A strong decrease in the fraction of mas- 
sive vs. low-mass stars is seen in Fig. 10. 
While 44 % of all stars with M > 1 M© 
in the cluster core, 7 < i? < 10 , have 
masses above 4 Mq (and M < 20 Mq), this 
fraction drops to 20 % already at a radius 
of 20 . At a radius of 33 , the contribu- 
tion of massive stars is below 5 %. Be- 
yond this radius, field contamination at the 
bright end skews the ratios again towards 
higher values. The steep decrease refiects 
the high-density concentration of massive 
stars in the cluster center, while at larger 
radii, i? > 30 , the cluster is dominated by 
the extended low-mass population. 

An even steeper drop is observed in the 
mass fraction of high- vs. low-mass stars 
(right panel in Fig. 10). While the total 



mass in high-mass stars exceeds the low- 
mass contribution by a factor of 1.8 close to 
the cluster center, meaning that high-mass 
stars contribute 64 % to the total mass, the 
ratio of the mass in high- vs. low-mass 
stars drops to 0.2 at i? = 33 , such that 
high-mass stars contribute only ~ 17% of 
the total mass at large radii while dominat- 
ing the mass in the cluster core. Note that 
in this estimate, the compact, massive core 
of the starburst is excluded due to satura- 
tion. The fiat core MF (F = —0.5) derived 
by Sung & Bessel (2004) for radii R < 6 
indicates a strong bias to high-mass stars 
in the core, such that including the massive 
core population would enhance the per- 
centage of mass contributed by high-mass 
stars in the center of NGC3603YC. 

As can be seen in Fig. 10, these vari- 
ations in the number and mass fractions 
with radius are a very efiicient tool to 
quantify mass segregation effects. 

8. From Orion to R136 - a struc- 
tural comparison 

Starburst clusters as massive and com- 
pact as NGC3603YC are rare in our close 
neighborhood. In the Milky Way, only the 
Arches cluster close to the Galactic Cen- 
ter is known to have comparable mass and 
density (Stohe 2003). Clark et al. 2005 ar- 
gue that the recently unveiled Westerlund 
1 cluster is comparable to Arches in core 
density and might have a total mass of 
10^ Mq under the assumption of a normal 
IMF, but only the most massive stars down 
to type BO are analysed so far, while the 
low-mass population is as yet unknown. 
In the literature, NGC3603 is extensively 
compared to the massive 30 Doradus star- 
forming complex in the Large Magellanic 
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Cloud (LMC), hosting the dense cluster 
surrounding the compact, star-like source 
R 136, which features a core density de- 
rived from O and WR stars comparable 
to the central NGC3603 starburst cluster 
(e.g., Moffat 1994). With a total mass of 
2 • 10^ Mq, the central cluster around R 136 
in 30 Dor is at the upper mass end of local 
young star clusters, and at the same time 
at the low-mass end of Milky Way globu- 
lar clusters. The large distance of 6 and 8 
kpc to NGC3603YC and Arches, respec- 
tively, allow only limited study of their low- 
mass population. This limitation applies 
to R136 and its surrounding cluster at 50 
kpc even more severely. Because of its 
proximity to the Sun, the Orion molecular 
cloud, on the other hand, is one of the best 
studied star-forming regions in the Milky 
Way. At a distance of only 450 pc, the 
Orion nebula is sufficiently close to resolve 
the embedded population down to the sub- 
stellar regime (Muench et al. 2002). The 
extended molecular cloud hosts the young 
Orion Nebula Cluster (ONC) containing 
stars with ages between 0.3 and 1 Myr 
in its center. The Trapezium system, a 
dense cluster of massive stars, forms the 
core of the ONC. It is particularly interest- 
ing for the comparison with NGC 3603 and 
Arches that the ONC provided a spatially 
well-resolved testcase where mass segre- 
gation could be analysed quantitatively 
(Hillenbrand 1997). The Trapezium clus- 
ter is with a core density and mass of 
about one order of magnitude lower than 
NGC 3603 YC at the low-mass end of mas- 
sive star-forming regions in the Milky Way 
(see Tab. 4), albeit an example of the mode 
of star formation we understand as normal 
in our local environment. As these star- 
forming regions play a major role in our 



understanding of the stellar IMF, we wish 
to compare their characteristics in the con- 
text of the influence of the star-forming 
environment on cluster formation. Tab. 4 
summarizes the known properties of these 
massive star-forming regions and their cen- 
tral clusters. Properties of the Antennae 
starburst clusters and Milky Way globular 
clusters are included for reference. 

8.1. Comparison of cluster charac- 
teristics 

8.1.1. NGC 3603 

The extent of NGC 3603 YC is esti- 
mated to a radius of 4.4 pc by Niirnberger 
& Petr-Gotzens (2002), beyond which it 
becomes indistinguishable from the stel- 
lar field population. The core radius is 
estimated to 0.2 pc from a King-profile 
fit to the density profile (Grebel et al., 
private communication). Within a ra- 
dius of ~ 1 pc, a total mass of 7000 M© 
is derived via extrapolation of the ob- 
served MF down to 0.1 Mq from a com- 
bination of IIST/WFPC2 core data and 
the ISAAC observations presented above 
(Stolte 2003). This estimate of the total 
cluster mass is a lower limit as it does 
not correct for the low-mass halo, which is 
difficult to quantify without independent 
field data. In the inner 0.12 pc {d = 6kpc) 
of the dense cluster core, six 03 stars and 
three WN6 components are identified by 
Drissen et al. 1995. 

8.1.2. R136 and its surrounding cluster 

As a consequence of the similarly com- 
pact central cluster and stellar mass con- 
tent in the inner 1 pc, NGC 3603 was sug- 
gested to be "a clone of R136" by Mof- 
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fat et al. 1994. Both clusters are with an 
age of 1-2 Myr equally young, and contain 
a massive 0-star population in their cen- 
ters. With 2 X 10^ M© (Hunter et al. 1995), 
the total mass in stars in the cluster 
around R 136 is a factor of 2-3 larger than 
in NGC3603YC. With the given stellar 
mass content, the central 30 Dor cluster 
was very productive in forming high-mass 
stars. Within the inner 0.2 pc, six 03 stars 
are classified by Massey & Hunter (1998), 
comparable to the six 03 stars confined 
to the core of NGC3603YC (r < 0.12 
pc). Despite this similarity in the imme- 
diate core, Massey & Hunter (1998) clas- 
sify 65 stars spectroscopically with types 
earlier than BO out to radii of 5.7 pc, of 
which 40 stars are of type 03, 4 are WN6 
stars, 6 are in the 03/WN6 transition 
phase to Wolf-Rayet stars, and one 03- 
03 binary is found, adding up to a total 
of 52 early O stars. In the 65 brightest 
stars with types BO and earlier spectro- 
scopically analyzed by Massey & Hunter, 
a strong bias to 03 stars is observed. A 
comparably large number of 03 stars is 
not known in any other local star-forming 
environment. Out of these, 21 03 ob- 
jects reside within 1 pc from R136, while 
no more 03/WN6 sources except for the 
nine very early type stars in the inner- 
most core of NGC3603YC are detected 
in the central 1 pc covered by the ISAAC 
data and further out to 4.2 pc (Moffat 
et al. 1994), the suggested cluster ex- 
tent (Niirnberger & Petr-Gotzens 2002). 
In NGC3603 YC, 84% (31) of the early O- 
star population {My > -5.0, M ~ 25 M©) 
is confined to the innermost 0.5 pc, and all 
0-stars are found within 1.15 pc from the 
cluster center. A comparable number of 
33 O and WN stars is classified by Massey 



& Hunter (1998) in the inner 1.1 pc of the 
cluster surrounding R136. However, this 
comprises only 63% of the 52 most lumi- 
nous stars in the inner 30 Dor cluster down 
to the same My limit, while the remainder 
is spread out to radii as large as 4.2 pc (see 
Moffat et al. 1994, Fig. 9). 

Comparable to the morphology ob- 
served in NGC 3603 YC, the cluster around 
R 136 displays a distinct dense cluster 
core surrounded by a stellar halo. The 
main difference between NGC 3603 YC and 
the core of 30 Dor is given by the spa- 
tial extent of the halo and its content 
of massive stars. While the dense clus- 
ter center around R136 is with a radius 
of 4.7 pc (Hunter et al. 1995) compara- 
ble to the total cluster radius of 4.4 pc 
derived for NGC 3603 YC by Niirnberger 
& Petr-Gotzens (2002), a halo of young, 
massive stars with Wolf-Rayet character- 
istics is detectable out to a distance as 
far as 130 pc from R136, beyond which 
it merges smoothly into the field popula- 
tion (Moffat et al. 1987). From the fact 
that only two more WNL stars are found 
in the rest of the LMC, Moffat et al. 1987 
argue that the observed young population 
of WN6/7 stars are members of 30 Dor. 
If R 136 and 30 Dor were formed in a sin- 
gle starburst, this renders the halo around 
R 136 a factor of 30 larger than the radial 
extent of NGC 3603 YC. 

8.1.3. Trapezium 

The ONC contains ~ 2000 M© within 
a radius of 2.5 pc (Hillenbrand 1997). 
The ONC appears as an extended halo 
of low- mass stars around the dense cen- 
tral Trapezium cluster. The core radius 
in the cluster derived from King model 
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fitting is 0.2 pc (Hillenbrand & Hart- 
mann 1998), identical to NGC3603YC. 
The central stellar number density inside 
the core is also very comparable, with 
Pc = 1.7 • 10"^ stars/pc^ in the Trapez- 
ium core and pc = 1.4 • lO^stars/pc^ 
in NGC3603YC. There is, however, one 
significant difference between both clus- 
ters - the central mass density is with 
~ IO^Mq/pc^ (Moffat et al. 1994) in 
the starburst about one order of magni- 
tude higher than the 2 ■ 10*^ Mq/pc^ in the 
Trapezium core. Given the large number of 
at least 40 0-stars in NGC3603YC (Dris- 
sen et al. 1995), versus a mere ~ 10 OB 
stars in the ONC, the difference in mass 
density is not surprising. NGC3603YC 
contains the entire mass found in the ONC 
only within its core. 

8.1.4- Arches 

The Galactic Center Arches cluster falls 
with a total mass of ~ 10^ M0 and cen- 
tral density of 3 • 10^ MqPc"^ inbetween 
NGC3603YC and the cluster surround- 
ing R136. As such, it is the most mas- 
sive, young compact cluster in the Milky 
Way for which a detailed MF analysis is 
available. The core radius is estimated 
to be < 0.24 pc from the half-mass ra- 
dius (Stolte 2003) comparable to both 
NGC3603YC and Trapezium, although 
the density profile cannot be used to de- 
rive a fitted core radius because it is heav- 
ily distorted by tidal disruption in the 
Galactic Center gravitational field. Due 
to the high stellar field density in the 
GC environment, the true extent and the 
population surrounding the cluster be- 
yond a radius of 1 pc is currently not 
known. A part of the cluster popula- 



tion has likely formed an extended tail or 
halo, as low- and intermediate-mass mem- 
bers are stripped off the cluster rapidly in 
the GC tidal field. Nevertheless, a com- 
pact, highly mass segregated high-mass 
core is observed in the center of a more 
extended intermediate-mass population, 
in structure comparable to NGC3603YC, 
but with a larger characteristic mass scale 
(Stohe et al. 2005; Stohe 2003). 

8.2. MF slopes in four clusters 

The central 4.7 pc of R 136/30 Dor 
exhibit a normal mass function with a 
slope of r = —1.3 ± 0.1 in the mass 
range 15 < M < 120 M©, identical to 
the slope of a Salpeter IMF. The mas- 
sive end, 30 < M < 120 Mq, of this MF 
is obtained from spectroscopy of the high- 
mass population (Massey & Hunter 1998). 
In the MF derivation, the 120 Mq upper 
mass limit is imposed from the truncation 
of stellar evolution models. From luminos- 
ity and - where available - spectroscopy of 
the brightest, highest-mass stars, Massey 
& Hunter suggest initial masses of up to 
150 Mq. Given the total mass of the clus- 
ter, these authors argue that this is in 
agreement with an untruncated normal 
mass function up to the highest masses. 

Hunter et al. (1996) study the radial 
variation in the MF slope in the mass 
range 2.8 < M < 15 Mq in four annuli 
with 0.11 < i? < 1.1 pc using HST/PC 
data. From inner to outer annuli, val- 
ues of -1.0 ± 0.4, -0.7 ± 0.4, -1.1 ± 0.3 
and —1.05 ±0.2^, respectively, are derived, 
showing no variation with radius within 



■''cf. Massey & Hunter (1998) for corrections of the 
values given in Hunter et al. (1996) 
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the uncertainties. The same behaviour 
was derived in the case of NGC3603YC 
in Sec. 6. The HST/PC data analysed by 
Hunter et al. 1995 only resolve the clus- 
ter beyond 0.1 pc. Mass segregation in the 
form of a flattened MF slope is observed 
in NGC3603YC only inside the core ra- 
dius of 0.2 pc (Sung & Bessell 2004), while 
segregation at larger radii is evidenced in 
the decrease in the mass fraction in high- 
mass stars (Sec. 7.2) and the truncation of 
the MF at the high-mass end (Sec. 6), but 
not in the slope of the MF. A mass seg- 
regated core biased to high-mass stars can 
thus not be entirely ruled out for R136, 
although the HST/PC image presented in 
Massey & Hunter 1998 (see their Fig. 2) 
suggestes that the center of R136 is well 
resolved with HST, and the MFs in Hunter 
et al. 1996 do not show a depletion at the 
high-mass end as seen in NGC3603YC. 
The mass function slope of F = —1.0 to 
-1.1 (2.8 < M < 15 Mq) observed in ap- 
proximately the same radial range, 0.1 < 
i? < 1.1 pc, agrees very well with the slope 
of F = -0.9 ± 0.1 (0.4 < M < 2OM0 in 
0.2 < i? < 1.0 pc) in NGC3603YC, both 
subject to a similar saturation limitation 
at the massive end. 

The mass function in the Trapezium 
system was investigated most recently by 
Muench et al. (2002). A slope of F = 
-1.2 ± 0.2 is found for 0.6 < M < 5 Mq, 
consistent with F ~ —1.3 found by Hil- 
lenbrand (1997) over a mass range 0.25 < 
M < 12 Mq, the upper hmit of the ONC 
IMF power law tail. These slopes are com- 
parable to the overall slope in R 136 and its 
surrounding cluster (2.8 < M < 120 Mq, 
Massey & Hunter 1998) despite the differ- 
ent mass ranges covered, supporting the 



idea of a universal power-law IMF above 
the turn-over mass. Again, no indication 
of a radially increasing slope is found out- 
side the cluster core (Muench et al. 2002). 

In the central parsec of the Arches clus- 
ter, the same slope of F = —0.9 ± 0.15 
as observed in NGC3603YC is obtained 
for the intermediate- to high- mass regime, 
4 < M < 65 M0. The low-mass regime 
of the cluster is not yet resolved. Al- 
though all of the above slopes barely agree 
within the errors, the low fltting uncer- 
tainties probably underestimate the true 
slope uncertainty, as photometric age or 
distance uncertainties are not taken into 
account. Thus, there may be a weak indi- 
cation of a flattened mass distribution in 
NGC3603YC and Arches, while the MFs 
in the clusters surrounding R136 and the 
Trapezium follow a standard Salpeter law. 

Although the spatial distribution of 
high-mass stars around R136 differs dras- 
tically from the morphology observed in 
Trapezium, NGC 3603 YC and Arches, the 
slope of the mass function is almost iden- 
tical. As already discussed in context of 
the radial variation and mass segregation 
in NGC 3603 YC, this supports the conclu- 
sion that the slope of the MF alone is not 
sufficient to understand mass segregation. 

Given the diversity of the star-forming 
environments considered, the slopes of the 
present-day MFs are in remarkable agree- 
ment, in accordance with a universal IMF 
slope. 

8.3. Primordial vs. dynamical seg- 
regation? 

In all four clusters, a heavily mass segre- 
gated core is observed. Mass segregation in 
young clusters was shown for the ONC and 
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Trapezium in Hillenbrand 1997 (see also 
Hillenbrand & Hartman 1998). Hunter 
et al. 1995 conclude from dynamical con- 
siderations that the core inside R 136/30 
Dor is entirely mass segregated for radii 
r < 0.4 pc, and their IMF suggests that 
at least 40 % of the total mass are con- 
fined within this radius. We estimate a 
half-mass radius of 0.24 pc for Arches, 
thus about 50 % of the mass are con- 
fined within the core, similarly indicated 
for NGC3603YC from the HST analy- 
sis (Stohe 2003; Sung & Bessell 2004). Do 
these heavily high-mass biased, mass segre- 
gated cores imply a preferential formation 
of high-mass stars in the cloud cores, and 
thus primordial segregation? 

Although we are not in the position 
to give a final answer to this question, 
there are several indications we can de- 
rive from the comparison of the known 
clusters. In all four clusters, the present- 
day MF rapidly approaches a slope close 
to a standard IMF outside the core ra- 
dius. In NGC3603YC we have seen that 
segregation mostly affects the high-mass 
stars, while the shape of the low-mass dis- 
tribution appears constant with radius. 
A similar behaviour is observed in the 
form of a constant MF slope in the clus- 
ter around R136 and Trapezium. In the 
Arches cluster the case is not clear due 
to the severe crowding and field contami- 
nation at masses below ~ 4Mq. Despite 
these similarities in the intermediate-mass 
MF, NGC3603YC and 30 Dor show sig- 
nificant structural differences in the spa- 
tial distribution of high-mass stars. While 
both NGC 3603 YC and Trapezium exhibit 
a halo of low-mass objects, the cluster sur- 
rounding R 136 is the only local starburst 



cluster that exhibits an extended halo of 
high-mass stars. Hunter et al. (1995) ar- 
gue that the dynamical relaxation process 
is completed around R136 in the inner- 
most 0.4 pc, with the relaxation timescale 
for high-mass stars being much smaller 
than the cluster age within this radius. 
This situation resembles the derivations 
for NGC 3603 YC (Stolte 2003) and Orion 
(Hillenbrand & Hartmann 1998) in the re- 
spect that a short dynamical relaxation 
time prohibits one to distinguish primor- 
dial and dynamical segregation in the clus- 
ter center. In the extended halo sur- 
rounding R136, however, the relaxation 
timescale is much larger due to the lower 
density and larger distance from the cluster 
center, such that stars formed in the outer 
regions of the cluster are not significantly 
segregated, but should still be found close 
to their birth positions. This implies that 
even the formation of the highest mass 
stars is not confined to the densest clus- 
ter region, as observed in the two Milky 
Way starburst clusters. Nevertheless, even 
the Milky Way hosts OB associations with 
lower stellar density capable of forming 
numerous O stars. The most prominent 
example discovered so far is the Cyg 0B2 
complex, which harbours 48 O stars with 
one star classified as early as 03 (Comeron 
et al. 2002, Hanson 2003). The young stel- 
lar population is extended over an area 
of ~ 20 pc x20 pc and lacks a compact 
core as opposed to the 1 pc into which the 
high-mass core of NGC 3603 YC is con- 
fined. Furthermore, populations hidden 
by extinction such as Westerlund 1 with 
an estimated core density of ~ 10^ MqPc^^ 
are just recently unveiled by deep NIR ob- 
servations. Westerlund 1 hosts at least 53 
early-type stars and a total mass of 10^ Mq 
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was suggested for the cluster (Clark et 
al. 2005), about one order of magnitude 
higher than estimated in NGC3603YC. 
Despite the comparable core density, O- 
type stars are found out to radii of 3pc 
from the cluster center in Westerlund 1, 
and the spatial distribution of massive 
stars resembles the environment around 
R136 more closely than NGC3603YC. In 
which respects is the star-forming environ- 
ment in the LMC 30 Dor complex different 
from the NGC3603 and Arches starburst 
environments? 

Several differences are observed be- 
tween 30 Doradus and the NGC3603 
molecular cloud. First of all, the to- 
tal gas mass in the 30 Dor complex is 
with (1..6) • 10^ Mo (Cohen et al. 1988, 
depending on the CO/H2 conversion fac- 
tor) at least one to two orders of mag- 
nitude larger than the 4 • 10^ Mq deter- 
mined with the same method in NGC 3603 
(Grabelsky et al. 1988). The lower metal- 
licity in the LMC is frequently employed as 
the origin for differences in the young stel- 
lar population. However, the metallicity in 
30 Dor has recently been shown to be 2/3 
of the solar value (Peimbert 2003), such 
that significant changes in the outcome of 
the star-forming process are unlikely due 
to the marginally lower metallicity. An- 
other significant difference was, however, 
suggested early: Werner et al. 1978 esti- 
mated an optical path length of 30 pc in 
30 Dor vs. 0.1 pc in Orion from the much 
lower far-infrared optical depth, indicating 
a low density of dust and molecular mate- 
rial in the giant cloud. If photons travel 
longer distances, the heating of molecu- 
lar cores may be much more efficient out 
to larger distances from the dense cloud 



center, where the first high-mass stars are 
likely to have formed. The suggested op- 
tical path length is on the order of 100 
times larger than in the Orion cloud. Tak- 
ing into account the quadratically decreas- 
ing radiation field intensity with distance, 
one expects the spatial extent of clusters 
in Orion to be about one order of mag- 
nitude smaller than in 30 Dor, consistent 
with the 2.5 pc radius of the ONC. If the 
gas/dust distribution in NGC 3603 is sim- 
ilar to Orion, this could explain the ~ 30 
times more extended halo of 130 pc around 
R 136 vs. 4.4 pc in NGC 3603. Observa- 
tions of core temperatures in the 30 Dor 
complex as well as a comparative study of 
optical path lengths in Milky Way GMCs 
would clearly be very valuable to probe 
this suggestion. 

A higher average cloud temperature is 
actually observed in the Galactic Center 
environment. Typical values of 70 K with 
a range of 30-200 K are measured in the 
GC (Morris & Serabyn 1996) vs. ~ 20 K 
observed in the cores of moderate star- 
forming environments. The intense radi- 
ation field in the GC may have a similar 
effect on core heating as the longer optical 
path length suggested for 30 Dor. As the 
path length in the GC is likely very small 
due to high gas and dust densities, the 
less efficient core cooling may locally cause 
the formation of cores with higher mass 
rather than extended formation of high- 
mass cores as in 30 Dor. A low cooling 
efficiency and intensive external heating 
may result in a higher core temperature 
when the opacity limit is reached and frag- 
mentation ceases, implying a lower core 
density and thus higher mass needed for 
gravitational collapse. A lower critical 
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density at which fragmentation stops re- 
sults in less fragmentation and a higher 
average core and final stellar mass ac- 
cording to recent numerical simulations 
(Jappsen et al. 2005). The present-day 
MF in the Arches cluster displays indi- 
cations for a high-mass turn-over around 
6M0 followed by a possible depletion to- 
wards lower masses (Stolte et al. 2005). In 
contrast to 30 Dor, where high-mass stars 
form at large distances from the cluster 
center, but the overall IMF appears nor- 
mal, an enhanced core temperature in the 
GC may have caused a locally top-heavy 
IMF in the Arches core. 

The Arches cluster is known to be 
tidally disrupted rapidly by the strong GC 
tidal field. Thus, for this cluster a skewed 
present-day MF may well have resulted 
from dynamical segregation alone. In 30 
Dor the spatial distribution of massive 
stars is unlikely to be caused by dynami- 
cal evolution, as the dynamical timescales 
in the outskirts of the cluster (beyond 
r = 1 pc) are larger than the cluster age 
(Hunter et al. 1995). Although NGC3603 
as a counterpart to the Arches cluster 
evolving in a moderate star-forming en- 
vironment is estimated to be dynamically 
significantly mass segregated (see Sec. 7), 
and will likely not survive on globular clus- 
ter timescales due to its compactness, the 
dynamical timescales in NGC3603YC are 
still at least one order of magnitude longer 
than in the GC environment (Stolte 2003). 
Thus, the differences between Arches and 
NGC 3603 YC can be caused by dynamical 
evolution, while the structural discrepan- 
cies between NGC 3603 and the central 30 
Dor cluster likely originate from intrinsic 
differences in the parental clouds. 



Hillenbrand & Hartmann (1998) per- 
formed a detailed dynamical analysis of the 
ONC and its core, the Trapezium, and the 
mass segregation studied by Hillenbrand 
(1997). The radial variation in the mass 
distribution observed in NGC 3603 YC and 
the ONC is comparable despite the fact 
that the starburst cluster hosts at least 
40 0-stars as opposed to only 3 0-stars 
in the Trapezium system. Hillenbrand & 
Hartmann derive cumulative mass distri- 
butions in different annuli for several mass 
bins. These CFs show a different shape 
for massive stars with M > 5 Mq, while 
the three mass bins between 0.3 and 5 Mq 
all show a remarkably identical behaviour. 
In the cluster core {R < 0.5 pc), how- 
ever, all four mass intervals follow each 
other closely, indicating strong segregation 
in the core, producing a bias to massive 
stars with respect to the low-mass popu- 
lation. The similarity over the wide range 
in masses suggests a fiat MF exclusively 
in the core of the ONC. At larger radii, 
mass segregation in the ONC is observed 
down to 1 — 2 Mq, but not below these 
masses. This is surprisingly similar to 
NGC 3603 YC, where the CFs of the PMS 
population up to masses of ~ 3 Mq are 
indistinguishable in different annuli, while 
for larger masses strong changes are ob- 
served (see Fig. 9). 

From the central stellar density, the 
number of stars in the core, the mean den- 
sity and the measured velocity dispersion 
of 2.5 km/s, Hillenbrand & Hartmann de- 
rive a relaxation time of 6.5 Myr in the 
ONC. This timescale is longer than the 
maximum age of ~ 1 Myr, suggesting 
that some amount of primordial segrega- 
tion had to be present. The same authors 
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also note that if only the massive core with 
~ 400 Mq is taken into account, the higher 
mean mass and stellar density causes a 
decrease to tr-eiax ~ 0.6 Myr, comparable 
to the age of the young cluster popula- 
tion. This means that in the Trapezium, 
too, observational evidence for primordial 
segregation is weak, as the massive stars 
may have migrated to the cluster center 
within the lifetime of the cluster. Bonnell 
& Davies (1998) perform detailed N-body 
simulations of ONC-type clusters, and find 
that the mean mass in the core as well as 
the ratio of high- to low-mass stars can 
only be recovered if the massive stars are 
initially placed close to the core, in agree- 
ment with primordial segregation. Given 
the high stellar mass and density in the 
core of NGC 3603 YC, some amount of pri- 
mordial segregation has likely contributed 
to shaping the cluster. 

A weak indication for primordial segre- 
gation is given by the low- mass MF. Both 
the MF as well as the CFs in NGC 3603 YC 
and Trapezium (Hillenbrand & Hartmann 199^ 
display increasing truncation at the high- 
mass end with increasing cluster center 
distance, but the shape of the low-mass 
distribution remains remarkably constant. 
From purely dynamical segregation a sys- 
tematic increase in low-mass stars with 
increasing radius and a corresponding de- 
crease in the apparent median (and char- 
acteristic) mass is predicted. While mea- 
suring the radially varying characteristic 
mass in the form of the turn-over in the 
MF in NGC 3603 YC has to await higher 
resolution observations, the mass fraction 
in low-mass stars increases with radius in 
NGC 3603 YC as shown in Fig. 10 con- 
sistent with expectations from dynamical 



evolution models. A detailed comparison 
with numerical simulations is needed to 
determine whether both the shape of the 
low- and high-mass MF as well as the in- 
creasing fraction of low-mass stars can be 
explained by dynamical segregation alone, 
or whether and how much primordial seg- 
regation is needed to reproduce the high- 
mass core along with the radially invariant 
properties observed in the low-mass halo. 
In summary, although we cannot finally 
distinguish between primordial vs. dy- 
namical segregation in the four clusters, 
we observed striking similarities between 
NGC 3603 YC and the ONC/Trapezium 
system, while structural differences are 
very prominent between NGC 3603 and the 
starburst surrounding R136 in 30 Dor as 
well as the Arches cluster in the GC envi- 
ronment. From these we conclude that the 
star- forming environment, despite surpris- 
ing similarities in the efficiency to produce 
the same relative fractions of high- and 
low-mass stars and thus the same IMF 
slope, influences the structural properties 
and thus the long-term survival of star- 
burst clusters severely, and therefore also 
the potential emergence of globular clus- 
ters from starburst environments. 

9. Summary and Conclusions 

From high-resolution near-infrared pho- 
tometry, the present-day MF in the Milky 
Way starburst cluster NGC 3603 YC is de- 
rived. Different treatments of the field star 
contamination, individual extinction cor- 
rection, and binary estimates are used to 
obtain the MF. We show that the slope of 
the MF is surprisingly robust against the 
derivation procedure. A slope of —0.91 ± 
0.15 is found in the mass range 0.4 < 
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M <20Mq for the integrated MF, in ex- 
cellent agreement with the optically de- 
rived MF by Sung & Bessell 2004, who 
obtain the same slope in the mass range 
1 < M < 100 Mq despite entirely different 
procedures employed to calculate the MF. 
This slope is identical to the slope found in 
the GC Arches starburst cluster, and only 
slightly flattened with respect to a Salpeter 
MF (r = -1.35). 

Radial variations are observed in the 
spatial distribution of stars, especially re- 
produced in the fraction of high- to low- 
mass stars and thus the characteristic mass 
at a given radius, and the depletion of mas- 
sive stars at the high-mass end of the MF. 
This radial variation is, however, not re- 
produced in the slope of the MF with val- 
ues of -0.9, -1.1, -0.8, -0.9, -0.9 in five equal- 
area annuli from immediately outside the 
cluster core out to radii of 2 pc. Although 
the slope is frequently employed as the only 
parameter comparing the MFs in different 
star-forming environments, it may not be 
the most significant value to use. The frac- 
tion of low-mass stars, the characteristic 
mass and the upper mass limit reflect mass 
segregation more closely. 

A detailed comparison between NGC 3603, 
30 Dor/R 136 in the LMC, Orion/Trapezium, 
and the Galactic Center Arches cluster 
suggests that a mass segregated core with 
an extended stellar halo may be a common 
cluster structure in a diversity of environ- 
ments. While the formation locus of the 
highest-mass stars appears confined to the 
cores of these two Milky Way starburst 
clusters, with the surrounding halo being 
comprised mainly of low- mass stars, stars 
as early as spectral type 03 seem capable 
of forming at large radial distances from 



the cluster center around R 136 in 30 Dor. 
Interestingly, the morphology in the Galac- 
tic star-forming complex Cyg 0B2 dis- 
plays similarities to the cluster surround- 
ing R136, and it would be interesting to 
compare the cloud temperature and spa- 
tial distribution of stars in both regions. 
This suggests that the star-forming envi- 
ronment, although producing comparable 
MF slopes, shapes the structural appear- 
ance of young star clusters, and thus in- 
fluences their long-term survival. While 
strong mass segregation is observed in the 
high-mass component of each cluster, and 
especially in the cluster core, the low-mass 
population appears constant in shape, de- 
creasing only as expected from the den- 
sity profile, over a large range of radii. 
This may indicate that some amount of 
primordial segregation had to be present 
in the clusters, even if dynamical segre- 
gation rapidly condensed the cluster core 
even more. Detailed N-body simulations 
adjusted to the derived cluster character- 
istiscs as were performed for the Trapez- 
ium cluster by Bonnell & Davies (1998) 
are required to distinguish primordial and 
dynamical segregation in NGC 3603 YC. 
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Fig. 1.— Ks 39 mill ISAAC image (3.'4 x 3.'4, North is up, East left). 

The selection of the main cluster center, i? < 33 (outer circle), as suggested by enhanced 
stellar density, and of the field population to the East of the cluster (white boxed area) are 
shown. Circles mark radial annuli with 7 , 20 , 27 , and 33 . 
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Fig. 2. — Density profile of NGC3603YC. Tlie dotted line indicates tlie limit where field 
star contamination exceeds 10%. 



Fig. 3. — Recovery fractions from artificial star tests. 

The recovery rates are shown for the cluster center and field including the polynomial fits 
used in the field star fraction and mass function incompleteness correction (left). The radial 
variation in the completeness due to varying stellar density is compared for the different 
annuli (right) for which the MF is calculated. Variations are significant only in the innermost 
annulus (7 < i? < 20 ). The huge loss of stars in the innermost 7 reflects the saturation 
and crowding problems in the cluster core, which is excluded from the MF derivation. 



Fig. 4. — Color- magnitude diagram of the inner 33 of NGC3603 YC. 

The unsubtracted cluster population is shown in the left panel, the observed field star popu- 
lation (not yet scaled to the same area) in the middle, and the subtracted cluster population 
used for mass function derivations in the right panel. Field stars were statistically subtracted 
in color-magnitude bins of 0.5 x 0.5 mag. The 50% completeness limit down to which mass 
functions are fitted is indicated as a dashed line. Note that the magnitude corresponding to 
50% completeness varies with radial distance from the cluster center. 



Fig. 5. — Selection of stars around the average MS (left) and PMS (right) isochrone. 
Stars with colors up to 0.25 mag blueward of the PMS and MS isochrones are included, and 
redward up to 0.35 mag in the case of MS stars and 0.5 mag in the case of PMS stars. The 
selected stars represent a typical sample used for the MF derivation (after field subtraction). 



Fig. 6.— Mass function of NGC3603YC. 

All MFs are derived from the combined MS and PMS population (Fig. 5) for an age of 1 

Myr, distance modulus of 13.9 mag and foreground extinction of Ay = 4.5 mag (PMS) and 

Av = 4: mag (MS). Dashed lines are individual fits for the PMS (0.4 < M/Mq < 4) and 

MS (4 < M/Mq < 20), and the solid red line is the combined MF fit. Top panel: MF 

calculated from simple star counts. Middle panel: MF calculated after field subtraction and 

individual dereddening. Bottom panel: Mass function derived including binary correction. 

Stars on the binary candidate sequence were shifted down by AJg = +0.75 mag, and 30% of 

the PMS stars selected randomly were treated in the same way. One equal-mass companion 

was included for each binary candidate. 

The similarity in the combined MF slope is striking despite the different procedures used. 
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Fig. 7. — Radial variation in the combined mass function. 

Three consecutive annuh with identical area coverage are chosen. In the combined MF a 
reasonable fit could be performed in the two inner annuli up to 16 Mq, while in the outermost 
annulus only the PMS contributes significantly to the stellar population, such that the fit 
proceeds only from 0.4 to 4Mq. Despite the observed decrease in the number of high-mass 
stars in the outer annuli, no clear steepening in the MF in the form of a steeper slope F 
reflects the depletion of high-mass stars at large radii. 



Fig. 8. — Cumulative mass distributions of NGC 3603 YC. 

The CF for the cluster selection 7 < i? < 33 is shown in the upper left panel. Theoretical 
mass distributions are overlaid (dashed lines) with slopes of —0.3 > F > —1.3 in steps of 
0.2, the latter corresponding to the Salpeter slope (solid line). All CFs are normalised to 
unity at M = 0.4 M©, and the absolute number of stars entering each CF is given in the 
upper right corner. The radial CFs are displayed in the subsequent panels, with annuli 
indicated in the upper right corner. In each radial CF, the cluster CF from the first panel 
is also shown for comparison. The two outer annuli extending to 65 indicate that the pre- 
main sequence population (M < 4:Mq) is not truncated even at large radii, while the main 
sequence population is strongly concentrated towards the cluster center. 



Fig. 9. — Radial cumulative mass distributions of NGC 3603 YC. 

This comparison between radial CFs of all annuli shows - in contrast to the binning-limited 
MF - a clear trend towards decreasing values of F with increasing radius in the inner three 
annuli. Beyond 33 , the increase at the high-mass end is most probably due to incomplete 
field subtraction. On the low-mass end, however, the mass distribution correponds very well 
to the cluster pre-main sequence distribution. Indeed, between all 5 annuli, no deviation in 
the distribution can be detected below 3 Mq except for a slight flattening in the cluster center 
(top curve). At masses above 3Mq, a strong decrease in the relative fraction of high-mass 
stars is observed, confirming the mass segregation observed in the MFs. 



Fig. 10. — Fraction of high- to low- mass stars vs. radius. 

The left panel shows the radial change in the number ratio of stars with M > 4 Mq to 
stars with 1 Mq < M < 4Mq (left). Large symbols represent the five radial annuli used to 
study MF variations, including individual incompleteness corrections for each annulus. The 
MS/PMS transition is chosen as mass limit between the high- and the low-mass population, 
as these two populations are known to be physically different, such that the ratio simultane- 
ously reveals the radial change in the number of MS vs. PMS stars. The strong decrease in 
the fraction of high-mass stars reflects the segregation of massive stars towards the cluster 
center. The right panel shows the radial variation in the mass fraction for the same criteria. 
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Table 1: Fraction of stars on binary candidate sequence. Only MS stars are selected. 

J g mag ^^binary ^^ single ^^ binary I y-'-^ single ' ^^ binary ) 

12 - 13 mag 9 19 0.32 

13 - 14 mag 12 38 0.24 
14 - 14.8 mag 15 42 0.27 

all 36 99 0.28 



Table 2: Mass Function derivations for NGC3603YC. 

Model number and description PMS + MS fit PMS fit MS fit 

no field subtraction, no dereddening r= -0.91 ±0.10 T = -0.78 ± 0.32 T = -0.75 ± 0.20 

with field subtraction, no deredding T = -0.87 ± 0.10 V = -0.72 ± 0.33 T = -0.72 ± 0.25 

with dereddening, no field subtraction T = -0.91 ± 0.10 T = -0.85 ± 0.30 T = -0.56 ± 0.20 

with dereddening and field subtraction r= -0.85 ±0.10 r = -0.73 ±0.32 r = -0.57 ±0.20 

binary candidates rejected r= -0.97 ±0.10 r = -0.73 ±0.33 T = -0.91 ± 0.12 

visible binary sequence corrected (MS and transition) r= -1.00 ±0.10 T == -0.73 ± 0.33 T = -0.95 ± 0.22 

visible binary seq. and statistical PMS correction 30% P = -0.91 ± 0.08 P = -0.67 ±0.26 r = -0.84 ±0.13 

2 Myr MS, all corrections P = -0.91 ± 0.08 P = -0.61 ± 0.21 P = -1.03 ±0.23 

3 Myr MS, all corrections (fit < 15 Mq) P= -0.84 ±0.09 P == -0.64 ± 0.22 P = -0.75 ±0.10 



Table 3: Radial variation of the mass function. 

annulus PMS + MS fit remark 



7" 


<R< 20" 


r = 


-0.89 ±0.14 


MS ± PMS fit 


20" 


<R< 27" 


r = 


-1.09 ±0.15 


MS ± PMS fit 


27" 


<R< 33" 


r = 


-0.76 ±0.21 


PMS fit only 


33" 


<R< 46" 


r = 


-0.99 ±0.26 


PMS fit only 


46" 


<R< 65" 


r = 


-0.85 ±0.23 


PMS fit only 
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Table 4: Comparison of young, massive clusters. 

cluster Mtotal extent r^ore Pcore 

Mp-, pc pc Mq pc"^ 



'0 



age 

Myr 



MF slope 

r 



ref 



Arches 


10^ 


1(?) 


0.2 


3-10^ 


2-3 


-0.9 ±0.15 


3,4 


NGC3603YC 


> 7-103 


4.4 


0.2 


10^ 


1-3 


-0.9 ±0.1 


4 


R136 


2-10^ 


4.7 


0.02 


5-10^ 


1-5 


-1.3±0.1 


1 


Orion 


10^ 


3 


0.2 


4-10^ 


0.3-1 


-1.2 ±0.1 


2 


Antennae 
















starbursts 


10^ - 10^ 


1-10 


? 


10^ 


1-20 


? 


5 


Milky Way GCs 


10^ - few 10^ 


few pc 


^ 1 


102 - 10^ 


10 Gyr 


-1.3 


various 



1 - Massey & Hunter (1998); 2 - Hillenbrand & Hartmann (1998); 3 - Figer et al. (1999); 4 - Stolte 
2003; 5 - Whitmore et al. (1999). Average Milky Way globular cluster parameters compiled from 



various sources. 
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